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Abstract. A new method for non-perturbative calculation of Green functions in 
quantum mechanics and quantum field theory is proposed. The method is based 
on an approximation of Schwinger-Dyson equation for the generating functional 
by exactly soluble equation in functional derivatives. Equations of the leading 
approximation and the first step are solved for ^-model. At d = 1 (anharmonic 
oscillator) the ground state energy is calculated. The renormalization program is 
performed for the field theory at d = 2, 3. At d = 4 the renormalization of the 
coupling involves a trivialization of the theory. 

1 Introduction and General Consideration 

During many years the construction of non-perturbative approximate solutions 
remains the urgent problem of quantum field theory. 

In present work a method for the construction of such approximations is pro- 
posed. The method is based on the approximation of the Schwinger-Dyson equa- 
tion for the generating functional by a simple equation in functional derivatives, 
which can be solved exactly. This solution is a foundation for the linear iterative 
scheme. Each step of the scheme consists in solving a closed system of integral 
equations. 

The leading approximation and the first step are investigated by the method 
for the <p£ model. At d = 1, when the model corresponds to the anharmonic 
oscillator, a formula for the ground state energy is obtained. At d = 2, 3 (super- 
renormalizable field theory) the renormalization of the leading approximation and 
the first step is performed. At d = 4 (strictly renormalizable case) the renor- 
malization of the coupling leads to a non-physical singularity of the amplitude. 
This is a reflection of the well-known triviality problem for the ^-theory in a 
non-perturbative region. 

Consider the theory of a scalar field (j)(x) in Euclidean space (x £ Ed) with the 
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action 

s(4>) = J dx{ l -(d^f + + \<p 4 }. (i) 

The generating functional of 2n-point Green functions can be written as 

oo 

G = £ G 2 V, (2) 

n=0 

where rj(x, y) is a bilocal source. The nth derivative of G over rj with the source 
being switched off is the 2n-point Green function G 2n . 

Using of the bilocal source is an essential point of the scheme presented, there- 
fore we shall consider the theory without spontaneous symmetry breaking with 
m 2 > 0. A brief discussion of the theory with a single source and with sponta- 
neous symmetry breaking is contained at the last section. At d > 2 corresponding 
counterterms should be included in the action for the cancellation of ultraviolet 
divergences. 

A simple way to obtain Schwinger-Dyson equation for the generating functional 
with bilocal source r/(x, y) is using the usual Schwinger-Dyson equation with single 
source j(x) at the presence of bilocal source r\ 

5j 6 (x) Sj(x) J 8j{x') 

Differentiating of equation (3) over j(y) and using the connection condition 

S 2 G _ 5G 
Sj(x)6j(y) 5r)(y,x) 

we obtain after switching off the single source j the following equation 

4A A - + (m 2 -d 2 )-——-2 / dx'r,(x,x') T ——-6(x-y)G = 0. 

drj{y,x)drj(x,x) drj{y,x) J drj(y,x') 

(4) 

which contains the source rj only. 

An idea of the presented iterative scheme is as follows: we shall consider "an 
equation with constant coefficients" as a leading approximation, i.e., equation (4) 
with the next-to-last term omitted. This term contains the source r? manifestly. 
The Green functions are the derivatives of G(rj) in zero and only the behaviour of 
G near rj = is essential, therefore such an approximation seems to be acceptable. 
The equation of the leading approximation will be 

^ $ 2 Gq ^ , 2 _ q2\ _ _ y)Go — (5) 

Srj(y,x)Sr](x,x) 5r)(y,x) 
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The term omitted contains the source and should be treated as a perturbation. 
Hence, the iteration procedure for the generating functional 



G = Go + Gi + • • • + G n + • • • (6) 
consists in the step-to-step solution of the equations 

4A / 6 \f ( ^ +(m 2 -d^-^--5(x-y)G n = 2 f dx>r,(x,x>)-^-. (7) 
5r]{y,x)Sr]{x,x) Sr]{y,x) J Sr)(y,x') 

The solution of the leading approximation equation (5) is a functional 

G = exp J dxdyr](y, x) A (x-y), (8) 

where Ao is a solution of the equation 

4A A (0) A (x-y) + (m 2 - d 2 ) A (x - y) = 5(x - y). (9) 

At d > 2 the quantity Ao(0) must be considered as some regularization. 

Equation (9) looks as the self-consistency equation, but differs in the coefficient 
at A: in the self-consistency equation the coefficient is three times greater. In this 
sense equation (9) is more similar to the equation for the propagator in the leading 
approximation of 1/iV-expansion. Certainly the similarity is completely superficial, 
since the principle of the construction of the approximation scheme is different. 

The solution of equation (9) is a free propagator 

Ao = ^ (10) 

with the renormalized mass /U 2 = m 2 + 4A Ao (0). The quantity Ao(0) is defined 
from the self-consistency condition. 

The propagator is a first derivative of G(rj) over the source r] : A = ^ |??=o- 
As can be easily seen, it is simply Ao for the leading approximation. 

Notice, that all higher Green functions of the leading approximation starting 
with four-point function G 4 = j% |^=o do not possess the correct connected struc- 
ture and, correspondingly, the complete bose-symmetry. The correct connected 
structure and other consequences of bose-symmetry (e.g., crossing etc.) will be 
restored in consecutive order at following steps of the iteration scheme. This is 
easy to see, for example, analysing the iteration scheme equations at A — > 0. Such 
a peculiarity of the iteration scheme is originated by bilocal source and is not 
something exceptional : as is well-known, the similar phenomenon appears also in 
constructing l/TV-expansion in the bilocal source formalism. 
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The first step equation for the generating functional G\ reads 

(PGi | (m 2 d 2 ) ^ 1 

Srj(y,x)6r](x,x) 5rj(y,x) 



4X T1 7^7 , + (m 2 -d 2 )—±--5(x-y)G l =2 f dx' V (x,x') A (x'-y)G 

di](y,x)dr){x,x) Or)(y,x) J 

(11) 

A solution of equation (11) is looked for in the form G\ = Pi(i])Go, where 

Pi = \Ft] 2 + Airy. Taking into account the leading approximation, equation (11) 

gives a system of equations for F and Ai: 

= S(x - y') A (x' - y) + 6(x - x') A (y - y'), (12) 

(fi 2 -8 2 ) Ai (x - y) + 4A Ai (0) A (x - y) + 4AF =0. (13) 

Equation (12) is the linear integral equation for the function F in the momentum 
space. A solution of the equation is 

F^l^j = AoOr - y') A (*' -y) + A (x - x') A (y - y') 

-4A J dxid yi A (x - xi) A (y - x 1 )K(x 1 - yi) A (yi - x') A (yi - y'), (14) 

where the kernel K is a solution of the equation 

K(x -y) = 2<5(x - y) - 4A J dx'L(x - x')K(x' - y), (15) 

and L(x — y) = Aq(x — y) is a single loop. Equation (15) can be easily solved in 
the momentum space. Its solution is 

k{p) = — — . (i6) 

W l + 4AL(p) V j 

Notice, that first two terms in formula (14) for F are the missed connected struc- 
ture of the four-point function of leading approximation. Hence the connected 
structure of four-point function is restored yet at the first step of iterations. 

To solve equation (13) for Ai is also quite simple. Taking into account the 
formulae above the solution can be written as 

A 1 (x-y) = -J dx'dy' A (x - x')E r (x' - y') A (y' - y), (17) 
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where 

S r (x -y) = [4A Ax (0) + 8A A (0)}S(x - y) + S(x - y), (18) 
S(x - y) = -(4A) 2 A (x-y) J dx'L{x - x')K(x' - y). (19) 

The quantity Ai(0) is defined by the self-consistency condition. 

At A -» 0, as is easy to see, A = A + Ai = AP ert + 0(A 2 ), where AP ert is 
the propagator of the perturbation theory, i.e. at small A the first step propagator 
reproduces correctly the first term of the usual perturbation theory in the coupling. 

In the general case, the solution of equation (7) for the n-th step of the iteration 
scheme is 

G n = P n ( V )G 0} (20) 

where P n is a polynomial in r/ of a degree 2n. Therefore at the n-th step the 
computation of Green functions reduces to solving a system of 2n linear integral 
equations. 



2 Anharmonic Oscillator 

At d = 1 the model with the action (1) describes the quantum-mechanical anhar- 
monic oscillator. The parameter m 2 corresponds in the case to a frequency of a 
harmonic oscillator described by quadratic terms. At d = 1 ultraviolet divergences 
are absent, and quantities Ao(0) and Ai(0) are finite. Consequently, formulae (8) 
- (19) are applied directly for the computation of Green functions. 

Since at d = 1 Ao (0) = j-, the self-consistency condition becomes the 
equation for a renormalized mass (or, more exactly, for a " renormalized frequency" ) 
^. 

^ = m 2 + — . (21) 

Here \i = \J~y? . The equation has the unique positive solution at all positive m 2 
and A. 

To calculate a ground state energy E one can use the well-known formula (see, 
for example, M, 0]) 

^ = G 4 (0, 0,0,0). (22) 

In this formula G 4 is the four-point (or, two-particle) function: G 4 = -gj^r \rj=o- 
With the formulae above for Green functions of the leading approximation and 
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the first step, we obtain the following formula for the ground state energy of the 
anharmonic oscillator: 



dE 1 1 2A 2A 

1\ ~ V + aT^ (1 " ^(M + 2^) 2 ^' ( } 



Here M = 4/i 2 + Integrating the formula with a boundary condition E \\ = q - 
taken into account, one can calculate the ground state energy for all values of the 
coupling: < A < oo. 

At A — > : E = m(| + f^ + 0(A 2 )) — the perturbation theory is reproduced 
up to second order. 

At A -> oo : E = eoA 1 / 3 + 0(A -1 / 3 ), and e = 0.756. The coefficient e differs 
by 13% from the exact numerical one EQ Xact = 0.668 (see, for example, ||). 

At A/m 3 = 0.1 the result of the calculation with formula (23) differs from the 
exact numerical one || by 0.8%, and at A/m 3 = 1 differs by 6.3%. 

Consequently, the first step formula (23) approximates the ground state en- 
ergy for all values of A with the accuracy that varies smoothly from (at A — > 0) 
to 13% (at A — > oo). Comparing these results with the results of other approxi- 
mate methods, in particularly, with the method of variational perturbation theory 
(VPT) |l|, ^] or the (^-expansion method [f|, we can see that this method gives best 
results for the intermediated coupling region A ~ 0.1m 3 . Really, at A/m 3 = 0.1 
the accuracy of first step of the calculations is better than the accuracy of fifth 
step of VPT (see Table 6 of ref. 0) and five times as higher in comparison with 
the usual perturbation theory in the order C(A 5 ). However, for the asymptotic 
region of strong coupling the accuracy of the given calculations is not so good, and 
other methods such as VPT give more exact results at A — > oo. It seems be likely 
that a combination of the method proposed with the methods of VPT type would 
gives good result for all coupling values. 



3 Super-renormalizable theory (d=2 and d=3) 

At d > 2 the action (1) should be added by counterterms for the elimination 
of ultraviolet divergences. First consider super-renormalizable theory (d = 2 and 
d = 3). It is sufficient to add counterterms of mass renormalization ^-(f) 2 and wave 
function renormalization ^(d^) 2 in the case. The Schwinger-Dyson equation has 
the form of equation (4) with the substitution 

m 2 ^m 2 + <5m 2 , d 2 -> (1 + 5z)d 2 (24) 
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There is no need to add a counterterm of wave function renormalization for the 
leading approximation, and the equation of the leading approximation will be 

4A A + (5ml + m 2 - d 2 )—^- - 5(x - y)G = 0. (25) 

or]{y,x)dri{x,x) dr]{y,x) 

The l.h.s. of iteration scheme equation (7) is of the same form as the l.h.s. of 
equation (25) is. At n > 1 the counterterms 5m 2 n and 5z n should be considered 
as perturbations. Therefore corresponding terms should be added to the r.h.s. of 
equation (7). So, the first step equation will be 

iX wk*) + « + ™ 2 - 82) ^y - i( * - ^ - 

dx r)(x, x )— j- - om 1 — + dzid — -. (26) 

6r)(y,x') 07](y,x) or](y,x) 

The normalization condition on the physical renormalized mass fj? gives us a 
counterterm of the mass renormalization in the leading approximation 

5ml = ii 2 -m 2 -4A A (0). (27) 

This counterterm diverges logarithmically at d = 2 and linearly at d = 3. 

For the first step of the iteration equation (12) for F remains unchanged, and 
its solution is described by the same formulae (14)-(16). The equation for Ai 
changes in correspondence with equation (26). Its solution can be written in the 
same form (17), but now for S r one gets 

S r = 4A Ai (0) + 5m\ - 25ml - bz x d 2 + S, (28) 

where S is given by formula (19). The normalization conditions 

E r (V) = 0, S;(V) = (29) 

give us the counterterms of the first step 

Szt = -£'(-v 2 ), 

5m\ = 25ml - 4A A 1 (0) - t(-fi 2 ) - n 2 t'(-fi 2 ). (30) 
The renormalized mass operator is 

£ r (p 2 ) = t(p 2 ) - £(V) " (P 2 + ^ 2 )S'( V), (31) 
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where, in correspondence with (19), 



t(p 2 ) = - ( 4A) 2 / A 2 / ^ 2Z(g2) • ( 32 ) 
V ' K ' J (2ir) d /i 2 + (p - g) 2 1 + 4AL(<z 2 ) v ; 

The counterterm <5zi is finite at d = 2, 3. The counterterm <5m 2 diverges as that 
of the leading approximation does, namely, logarithmically at d = 2 and linearly 
at d = 3. As the simple loop L(p 2 ) behaves for p 2 — > oo as ^ log ^ at d = 2 
and —4= at d = 3, the integral (32) for £(p 2 ) converges at d = 2 and diverges 



logarithmically at d = 3. Surely the renormalized mass operator (31) is finite in 
any case. 

A part of "redundant" subtractions in formula (31) completely unclear from 
the point of view of the usual perturbation theory divergences. This part clears at 
the strong coupling limit in formulae (31)-(32). Really, at A — > oo 

tip 2 ) = const + ( -—(-t^ — r ^~^r +0(t)- (33) 

V1 ' J (2vr) d n 2 + (p - q) 2 L(q 2 ) V V ; 

The integral in equation (33) diverges quadratically at d = 2, 3. Hence, it becomes 
clear that the "redundant" subtractions conserve the finiteness of renormalized 
theory in the strong coupling limit. 

4 Strictly renormalizable theory (d=4) 

At d = 4 besides the renormalizations of the mass and the wave function a coupling 
renormalization is necessary. Therefore simultaneously with the substitution (24) 
the substitution A — > A + 5 A is also needed in the Schwinger-Dyson equation (4). 
The leading approximation equation will be 

Due to the presence of the counterterm 5X the normalization condition on the 
renormalized mass a 2 for the leading approximation becomes a connection between 
counterterms 5rriQ and 5Xo 

5ml + 4(A + <5A ) A (0) = a 2 - m 2 . (35) 

As we shall see below, counterterm 5\q (and, consequently, 5m 2 ,) will be fixed on 
the following step of the iteration scheme. 
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The first step equation will be of the form (26) with the substitution A — > 
A + 5\q in the l.h.s. and with an additional term — 4<5Ai • 5 2 Go/5r)(y, x)Sr/(x, x) in 
the r.h.s. The equation for F will differ from equation (12) by the substitution 
A — > A + 5\q only. Therefore the formulae for its solution also will differ from 
formulae (14)- (16) by the same substitution. At d = 4 the single- loop integral 
L(p 2 ) diverges logarithmically, and the renormalization of the coupling is necessary. 
Let define a two-particle amplitude — the amputated connected part of the four- 
point function 

A = A 1 A 1 F con A 1 A \ (36) 

Here a multiplication by Aq 1 is understood in the operator sense. F con is the 
connected part of F. It is easy to see that the amplitude depends from a variable 
V = Px + Py only and has the form 

A(p 2 ) = 8(A + <5A °l . (37) 

1 + 4(A + SX )L(p 2 ) 

Define a renormalized coupling A r as a value of the amplitude in a normalization 
point 

A(M 2 ) = -8A r = 8(A + S M . (38) 

V ' 1 + 4(A + <5A )L(M 2 ) V ; 

From equation (38) one obtains a counterterm of coupling renormalization 

<5A = -A + -~ (39) 

l-4A r L(M 2 ) V ' 

and renormalized amplitude 

A(p2) = 8 A (40) 

' l + 4A r L r (p 2 ;M 2 ) V ; 

where L r (p 2 ;M 2 ) = L(p 2 ) — L(M 2 ) is a renormalized loop that possesses a finite 
limit at the regularization removing. 

Taking the renormalization of the two-particle amplitude in such a manner, 
one can solve the equation for Ai and renormalize the mass operator in corre- 
spondence with the general principle of normalization on the physical mass (see 
equation (29)). But in four-dimensional case one gets an essential obstacle. At the 
regularization removing, <5Ao — > —A. It is evident from equation (39). Therefore 
the coefficient A + 5Xq in the leading approximation equation (34) vanishes. The 
same is true for all the subsequent iterations. The theory is trivialized. One can 
object that an expression 

^ °^ Srj(y,x)6r](x,x) ^ ^ 
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is really an indefinite quantity of • oo type, and the renormalization is, in the 
essence, a definition of the quantity. But it does not save a situation in this case 
since the renormalized amplitude (40) possesses a non-physical singularity in a 
deep-euclidean region (it is a well-known Landau pole). The unique noncontra- 
dictory possibility is a choice A r — ► at the regularization removed. This is the 
trivialization of the theory again. This trivialization appears almost inevitably in 
an investigation of ^-theory beyond the perturbation theory and is a practically 
rigorous result (see pi). Notice, that contrary to perturbation theory which is 
absolutely nonsensitive to the triviality of the theory, the method proposed leads 
to the trivialization already at the first step. 



The method considered above is based essentially on the bilocality of the source. 
Since the bilocal source produces 2n-point functions only, the method cannot be 
applied in its present form to a theory with spontaneous symmetry breaking when 



For description of the spontaneous symmetry breaking it is necessary to switch on 
a single source j(x), i.e. to consider Schwinger-Dyson equation (3). 

Consider a theory with the single source j and with the bilocal source switched 
off: T) = 0. Schwinger-Dyson equation for the generating functional G(j) is 



Let apply to equation (43) the same idea of approximation by the equation 
with "constant" coefficients, i.e. consider as a leading approximation the equation 



5 Single Source 



< | <p | >^ 0. 



(42) 




(43) 




(44) 



Then an iteration scheme will be described by the equation 




(45) 



The leading approximation equation (44) has a solution 




(46) 
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Surely v does not depend on x in a translational invariant theory. Therefore 
an equation for v will be 

4Au 3 + m 2 v = 0. (47) 

At m 2 > 0, A > the equation has unique real-valued solution v = 0, that 
corresponds to leading approximation Go = 1- Iteration scheme (45) based on this 
leading approximation coincides with the perturbation theory in the coupling - 
the leading approximation is too simple and does not contain any nonperturbative 
effects. 

At m 2 < besides of this solution the following real-valued ones exist 
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-*V-«' (48) 

which correspond to spontaneous breaking of a discrete symmetry (P-parity) of 4> - 
theory. A calculation of the ground state energy based on the formulae of equation 
(22) type demonstrates that the state with spontaneous symmetry breaking is 
energetically preferable and so it is a physical vacuum of the theory at m? < 0. 
There is the way to describe the leading non-perturbative effect, i.e. spontaneous 
symmetry breaking, by the method. 

A first step equation with counterterms will be 

a \ ,2 a2\ 

4A » + (m • 



5p(x) Sj(x) 

■I \n as\ ^ 3( ^0 s 2 ^0 . g fl 2 ^0 

j(x)Gq — 4oAi „ — om\ + bz\d 



5p(x) Sj(x) 5j(x) 

= (j(x) - 45X lV 3 - v5m\)GQ. (49) 

A solution of equation (49) should be looked for as G\ = P\(j)Go, where 

Pi{j) = \ ^i J 2 + The equation for Ai with taking into account formulae of 

leading approximation (46)- (48) has a solution 



A, = (50) 



where fi 2 = —2m 2 > 0. The reconstruction of the vacuum leads to the cor- 
responding reconstruction of the one-particle spectrum. All the picture exactly 
corresponds to a description in the effective potential language, but the concep- 
tion of the effective potential is not used at all. Subsequent iterations leads to 
the renormalized perturbation theory over the physical nonsymmetrical vacuum. 
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A remarkable feature of the scheme is an absence of symmetry breaking countert- 
erms even at intermediate steps of calculations. For the ultraviolet divergences 
removing the counterterms 5m 2 , SX and 5z are sufficient. 

In conclusion notice, though at m 2 > values of v in equation (48) are imagine, 
the corresponding real-valued solutions of equation (44) exist, for example 

Go = cosjii; J dxj(x)}, (51) 

where w 2 = m 2 /4A. At first step of the iteration such a leading approximation 
gives tachyons and so it is physically unacceptable. It is possible that similar 
solutions can be useful for an investigation of the problem of spontaneous sym- 
metry breaking in ^-theory with m 2 > (see, for example [||). Of course, the 
computational scheme should be modified in that case. 
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